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Abstract 

A detailed study is presented of elastic WW scattering in the scenario that there are no new 
particles discovered prior to the commissioning of the LHC. We work within the framework of the 
electroweak chiral lagrangian and two different unitarisation protocols are investigated. Signals 
and backgrounds are simulated to the final-state-particle level. A new technique for identifying 
the hadronically decaying W is developed, which is more generally applicable to massive particles 
which decay to jets where the separation of the jets is small. The effect of different assumptions 
about the underlying event is also studied. We conclude that the channel WW — » jj + Iv may 
contain scalar and/or vector resonances which could be measurable after 100 fb _1 of LHC data. 
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1 Introduction 



It is quite possible that no new particles will be discovered before the start of the Large 
Hadron Collider (LHC). Nevertheless, it is certain that new physics must reveal itself in 
or below the TeV region and it is likely that the LHC will be able to study this new physics 
in some detail. Precise data collected at LEP, SLC and the Tevatron interpreted within 
the Standard Model (or supersymmetric extensions of the Standard Model) suggest that 
this new physics should manifest itself as a higgs boson with mass less than around 200 
GeV However, such a limit is model dependent and it is possible for there to be no 
light scalar particle at all [||, §). 

The scattering of longitudinally polarised vector bosons via the process WlWl — > WxPFlq 
is particularly sensitive to the physics of electroweak symmetry breaking for it is in this 
channel that perturbative unitarity is violated at a centre-of-mass energy of 1.2 TeV. Thus 
we know that interesting physics must emerge before then. In the absence of a light higgs, 
or any other new physics, below some scale A, one can develop a quite general, model 
independent treatment of physics well below A. This treatment is underpinned by the 
electroweak chiral lagrangian (EWChL) In this paper we investigate sensitivity at 
the LHC to new physics within the EWChL. 

The process WW — > WW at high energy hadron colliders has been studied previously, 
usually in the context of searches for a heavy Higgs (for an overview see || || |7| |§). 
The ZZ decay modes constitute the principal discovery channel for Higgs masses above 
160 GeV or so, and the WW channels become important around 600 GeV. Within the 
chiral lagrangian, it has been usual to focus on leptonic decay modes of the gauge bosons 
in order to reduce hadronic backgrounds ||, [ji]]. In this paper we focus on the more 
complicated semi-leptonic final state. Cuts developed in previous studies || 0, [ll], [l2|, 13] 



are re-examined as a tool for measuring the cross-section differential in the WW invariant 
mass in the general case (i.e. with no assumption as to the presence or otherwise of a 
resonance). A novel technique for identifying the hadronic decays of boosted massive 
particles using the longitudinally invariant kr algorithm jl4| is introduced, and applied 
to identification of the hadronically decaying W . We also examine the sensitivity of 
the cuts and reconstruction methods to current simulations of the underlying hadronic 
activity. 

The paper is set out as follows: The EWChL formalism is introduced in Section 2, and in 
Section 3 we discuss the unitarisation of the scattering amplitude for WlWl — > WlWl- 
Unitarisation often leads to the prediction of resonances. We investigate the model 
dependence of such predictions and the nature of the resonances (scalar or vector) by 
looking at two different unitarisation protocols. In Section 4 we present parton level 
predictions for the WW production cross-section at the LHC for a variety of possible 
scenarios. The goal for the LHC will be to distinguish between these different physics 
scenarios. To study the potential for this, we have implemented the general formalism of 
the EWChL in the pythia Monte Carlo program (T^j. Sections 5 to 8 cover our analysis 

x We often use the symbol W to denote both W and Z bosons. 
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of both signal and background. We succeed in reducing the background to manageable 
levels using a variety of cuts which are discussed in detail. Section 9 contains a summary 
and conclusions. 

2 The Electroweak Chiral Lagrangian 

In the EWChL approach, new physics formally appears in the lagrangian via an infinite 
tower of non-renormalisable terms of progressively higher dimension. However, correc- 
tions to observables arising from the new physics can be computed systematically by 
truncating the tower at some finite order. This is equivalent to computing the observable 
to some fixed order in E/A where E is the relevant energy of the experiment. 

The breaking of electroweak gauge symmetry already informs us that the scale of this 
new physics should be around v = 246 GeV and the degree of symmetry breaking dic- 
tates that our lagrangian should involve three would-be Goldstone bosons (7?). Moreover, 
experiment has told us that after symmetry breaking there remains, to a good approx- 
imation, a residual global SU(2) symmetry (often called custodial symmetry) which is 
responsible for a /^-parameter of unity (p = M^-/(M^ cos 2 0\y)). I n chiral perturbation 
theory the residual SU (2) symmetry is the result of the breaking of a global chiral sym- 
metry, SU(2)l x SU{2)r. With these constraints, there is only one dimension-2 term 
that can be added to the standard electroweak lagrangian with massless vector bosons. 
It is 

£( 2 ) = ^(DpUDHUi) (1) 
where (• • •) indicates the SU(2) trace, and 

U = cxp (i^f} (2) 

(r are the Pauli matrices). This term contains no physics that we do not already know. It 
is responsible for giving the gauge bosons their mass (this is easiest to see in the unitary 
gauge where U = 1). 

At the next order in the chiral expansion, we must include all possible dimension-4 terms. 
There are only two such terms that will be of relevance to us. They are 

£ {4) = aidD^UD^)) 2 + a 5 ((D^UD^}) 2 (3) 

where CI4 and 05 parametrise our ignorance of the new physics and they are renormalised 
by one-loop corrections arising from the dimension-2 term. There are a number of addi- 
tional dimension-4 terms that can arise. However they generally contribute to anomalous 
trilinear couplings between vector bosons. In this paper we focus only on the quartic 
couplings. In the particular case of the Standard Model with a heavy higgs boson of 
mass mii, 0,5 = v 2 /(8mjj) and 04 = before renormalisation, whilst for the simplest 
technicolor models 04 = — 2as = Ntc/{^^ 2 )- 
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To date, other than fixing the scale v the main constraint on the parameters of the 
EWChL come from the precision data on the Z . Bagger, Falk & Swartz have shown 
that the EWChL can be accommodated without any fine tuning for A all the way up to 3 
TeV (general arguments based on unitarity indicate that A < 3 TeV) Q. They show that 
the Z° data constrain the couplings associated with a dimension-2 custodial symmetry 
violating term and a dimension-4 term which contributes to the electroweak parameter 
S. There are however no strong constraints on 04 and a§ and in this paper we assume 
that they can vary in the range [-0.01,0.01] ||i~6| . 

To one-loop, the EWChL yields the following key amplitude (p is the renormalisation 
scale) [§: 

1 10s 2 + 13(£ 2 + u 2 ) 



s 4 

A(S,t,U) = -n+^ 
1 



2a 5 {n)s z + a A {p)(t z + u z ) + 



(4vr) 2 72 



967r 2 v j 



t(s + 2t) log(-^) + u(s + 2u) log(-^) + 3s 2 log(-^ 



(4) 



■A 

in terms of which the individual WlWl — ► WlWl isospin amplitudes can be written: 

A (s,t,u) = 3A(s,t,u) + A(t,s,u) + A(u,t,s) (5) 
Ai(s,t,u) = A(t,s,u) - A(u,t,s) (6) 
A 2 (s,t,u) = A(t,s,u)+A(u,t,s). (7) 

Equation @ is derived assuming the Equivalence Theorem wherein the longitudinal W 
bosons are replaced by the Goldstone bosons |fT^| . This approximation is valid for energies 
sufficiently large compared to the W mass. 

In addition, @ is useful only for energies well below A, where the effects of the new 
physics manifest themselves as small perturbations. At the LHC, we will be hoping to 
see much more than small perturbations to existing physics. For example, we might see 
new particles associated with the physics of electroweak symmetry breaking. It would 
be very useful if we could in some way extend the domain of validity of the EWChL 
approach to at least address the physics that might emerge around the scale A. To a 
degree, this can be done by invoking some unitarisation protocol which ensures that @ 
develops a high energy behaviour that is consistent with partial wave unitarity [Q. In 
the next section, we will consider protocols that do not spoil the one-loop predictions of 
the EWChL at lower energies. Such an approach has met with some success in extending 
studies of chiral perturbation theory in QCD |^9[| . We will focus on two unitarisation 
protocols: the Pade protocol and the N/D protocol. 



3 Unitarisation 

The amplitude in the weak isospin basis, Aj, can be projected onto partial waves, t/j, 
with definite angular momentum J and weak isospin I: 

1 r 1 

tlJ = 64^ J 1 d(c ° S ^ Pj(c ° S ^ Al( - S ' *' U) (8) 
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where 9 is the centre-of-mass scattering angle. The WW scattering system can have 
I = 0, 1, 2 and Bose symmetry further implies that only even J are allowed for I = and 
2, while only odd J are allowed for I = 1. Subsequently we consider the three amplitudes 
ioo^lljfeo- The higher partial waves are strictly of order s 2 /v 4 but they are numerically 
small and we neglect them. 



Writing tjj = tfj + tjj H , the first two terms of the expansion are given by | 10|] : 



*(2) 
t 00 



167rt> 2 



(9) 



'oo 



647T1> 4 



16(lla 5 (/x) + 7a4(/i)) J_ / 101 - 501og(s//x 2 ) + ^ 



16tH 



f (2) 
'11 



967nr 



(10) 



fill 



£ (4) 



96vTO 4 



4(o4(/i)-2o 5 (/i)) + 



,(2) 



1 ( 1 i7T 
+ 



16vr 2 V9 6 



1 20 



32irv 2 



l 2() 



64tto 4 



Using 



32(o 5 (/i) + 2o4(m)) , 1/91 201og(s//* 2 ) . N 

n h ^vr 

3 16vr 2 I 18 9 



Ai(8,t,u) = 32vr J^(2J + l)t/jPj(cosl 
j=o 



we have (neglecting higher partial waves) 

A (s,t,u) = 32irt 00 
Ai(s,t,u) = 32-7T 3in cos # 
A 2 (s,t,u) = 327r£ 2 o- 

In terms of these amplitudes we can write 



A(W + W~ -» 

A(TU+1U~ -> ZZ) 
A(ZZ -> ZZ) 
A(WZ -» PFZ) 



1 



1 „ 1 , 
-A + -A 1 + -A 2 

3 2b 

1 A 1 A 

-An A 2 

3 u 3 

1 A 2 A 

-An + -A 2 

3 u 3 
1 , 1 , 

2* + 2^ 2 

An. 



(12) 
(13) 

(14) 
(15) 



(16) 



(17) 
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The differential WW cross-section is 



da \A(s,t)\' 



dcos9 32ttM^ w ' y ' 

To obtain the cross-section for pp — > WWjj + X we need to fold in the parton density 
functions, fi(x,Q 2 ), and the luminosity: 

da x ^ f 1 f 1 dzi cfo 2 2 w / fl/f 2 \ <*£ww f l da 



dM^ w JM^ w /sJM^ w /{x l s)Xxx 2 s pp J ' lx ~"^ w ' J: ' y ~"^ w/ dr J-idcos 



(19) 



where \/s pp is the centre-of-mass energy which we take to be 14 TeV, as appropriate for 
the LHC, 

dLwW ' f a ^ 1 [(l+r)ln(l/r)-2(l-r)] 



dr V 4-7T sin 2 Ow J T 



for incoming W bosons [20] and r = M^ w / (x±X2S 



VP) 



The Pade protocol Otherwise known as the Inverse Amplitude Method, this is a simple 



unitarisation procedure, and is widely employed |10], ^1], |2^|. Elastic unitarity demands 

that for s > the imaginary part of the amplitude is equal to the modulus squared of 
the amplitude, which implies 

tj} = Re(tj}) - i. (20) 
To the accuracy in which we work, we can write 

t (2) 

tIJ = TT^fen (21) 

which has the virtue that it satisfies the elastic unitarity condition identically. We stress 
that this method of unitarisation leads to an amplitude that is equivalent to the one-loop 
EWChL calculation modulo higher-orders in s/v 2 . 

Having unitarised the amplitude it is natural to ask what the consequences are. Typically, 
the partial waves develop resonances which serve to implement the demands of unitarity; 
this is the role played by the Higgs boson in the Standard Model. The position and nature 
of the resonances depends critically upon the unitarisation protocol and we investigate 
an alternative protocol in the following subsection. At high enough energy, the partial 
waves effectively lose all memory of the underlying chiral perturbation theory and their 
nature is driven solely by the choice of unitarisation protocol. We therefore rely on our 
unitarisation protocol to provide us with some feeling for the pattern of lowest lying 
resonances which may be observed at future colliders. 
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Figure 1: Map of the parameter space as determined by the Pade protocol [10|. The 
small triangle in the centre is the region of no resonances and the region below 
the dotted line is forbidden. Also shown are the points corresponding to the 
various scenarios considered in the text. 



Resonances are found whenever the corresponding phase shift passes through tt/2, i.e. 
when 

cot 5/ j = Re(ijj) = 0. 
In the Pade approach we can solve this equation to obtain the corresponding masses and 

4v 2 



widths 1 10]. For scalar resonances: 



m| 



^(llaefa) + 7o4(m)) + — 1 ™=™z£& M 



16vr 2 



(22) 



and 



For vector resonances: 



and 



mt 



my 



16ttv 2 

„,2 



4(a 4 (/x)-2a 5 (/i)) + T ^i 



(23) 



m 



V 



V 



96-7rf 2 

There are no resonances in the isotensor channel, i.e. from £20- There is however a region 
of parameter space where the phase shift passes through —ir/2. This would violate 
causality and so we are forced to forbid such regions of parameter space. It occurs when 



32 



1 



/273 20 s N 
f(«5(M) + 2a 4 ( M )) + I ^2( 1I - ¥ log^l<U. 



(24) 
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A map of the 04 — 05 parameter space showing the corresponding resonance structure is 
presented in Figure [[]. We fix fx = 1 TeV and, using equations ( |22| ) and (|23|), we define 
the regions to contain a resonance of the specified type with mass below 2 TeV. The blue 
points labelled TC and SM correspond to the naive Ntc = 3 technicolor (TC) and 1 
TeV Standard Model Higgs (SM) models. 



The N/D protocol This provides our alternative to the Pade protocol. This method 
ensures that the amplitude has improved analytic properties in addition to satisfying 
partial wave unitarity and matching the one-loop EWChL calculation. The right-hand 
cut is placed wholly into the denominator function, D, while the left-hand cut is encap- 
sulated in the numerator function, N, i.e. analyticity and unitarity demand the following 



relations fl22j, 23, U 



Im(t/j(s) _1 ) = -1 s > (25) 

ImL>(s) = s<0 (26) 

Im N(s) = Dim t u (s) s<0 (27) 

ImJV(s) = s>0 (28) 

«»W = £gf (29) 
Following Oiler, we define the following function to contain the right-hand cut at s = M 2 



where 



1. t 



7T 



g{s) = _i og (__l (so) 



where M is an unknown parameter. The N/D unitarised partial wave amplitude is then 
written 

tu(s) = 7~4#n ( 31 ) 

1 + g(s)X IJ (s) 

where 

Xjj(s) = tf)(s) + t^(s)+g( S )(tf)(s)f. (32) 



The amplitude thus defined has been constructed so as to satisfy (|25|) and (|28|) identically 
whilst ( |26| ) and (p7[ ) are satisfied to one-loop in chiral perturbation theory. Note that 
the contribution to ImD(s) for s < is beyond the one-loop approximation. 

In Figures Q to ^ we show curves of constant resonance mass, varying from 600 GeV to 
2 TeV in steps of 100 GeV, as a function of the appropriate combination of 04(1 TeV), 
05(1 TeV), and M. The horizontal lines obtained using the Pade protocol are tangent 
to the corresponding N/D contours. Over large regions of parameter space, the two 
protocols yield similar results. However, the N/D method predicts a larger region without 
resonances, indeed for M below around 1 TeV there are no resonances at all. Referring 
back to Figure [l], we see that as M increases the lines which define the scalar and vector 
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Scalar Channel 




Figure 2: Contours of constant scalar resonance mass in steps of 100 GeV. The horizontal 
lines are obtained using the Pade protocol and the curved lines are obtained 
using the N/D method which depends upon the parameter M. 



Vector Channel 




Figure 3: Contours of constant vector resonance mass in steps of 100 GeV. The horizontal 
lines are obtained using the Pade protocol and the curved lines are obtained 
using the N/D method which depends upon the parameter M. 
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Figure 4: Contours of constant isotensor resonance mass in steps of 100 GeV. The curved 
lines are obtained using the N/D method which depends upon the parameter 
M. The horizontal lines and dotted curves populate the unphysical region of 
parameter space, see the discussion in the text. 

regions move slowly outwards. From Figure [| we see that for M above ~ 1 TeV the region 
excluded in the Pade protocol is not excluded in the N/D protocol leading instead to a 
region without any resonances. For M below ~ 1 TeV, there is a region excluded by N/D 
unitarisation (for the same reason as in the earlier Pade case)|. The line delineating the 
forbidden region in Figure [l] thus moves slowly downwards as M decreases. Note that we 
do not know the natural value for M, e.g. it can be much smaller than 1 TeV. Finally, we 
note from Figure |] that the N/D method does allow for the existence of doubly charged 
resonances. 



4 Parton Level Predictions for the LHC 

In this section the parton level predictions for the process pp — > W + W~jj + X at 14 
TeV centre-of-mass energy are presented for the 5 different choices of 04 and 05 shown 
in Table |]. In the Pade approach these choices produce a 1 TeV scalar (scenario A), a 1 
TeV vector (scenario B), a 1.9 TeV vector (scenario C), a 800 GeV scalar and a 1.4 TeV 
vector (scenario D), and a scenario with no resonances (scenario E). The green points 
labelled A-E on Figure || correspond to the 5 scenarios we consider. Throughout this 



paper the CTEQ4L |25|] parton density functions as implemented in PDFLIB [26] are 



used, evaluated at the WW centre-of-mass energy (Mww)- The renormalisation scale 
2 The scalar and vector sectors have exclusion regions similar to the tensor sector. 
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Scenario 


a 4 (l TeV) 


o 5 (l TeV) 


A 


0.0 


0.003 


B 


0.002 


-0.003 


C 


0.002 


-0.001 


D 


0.008 





E 









Table 1: Parameters for the five scenarios which we consider. 

0.10 I i i i 1 i i i i i i i i i i i 1 r 




500 1000 1500 2000 

M„ (GeV) 

Figure 5: Parton level cross-section for Scenario A. We compare the Pade result (solid 
line) with the N/D results for M = 10 3 GeV(dashed line), M = 10 4 GeV 
(dashed-dotted line) and M = 10 5 GeV (dotted line). 

is fixed to 1 TeV. The differential cross-section do / dMy/w for each of scenarios A-E are 
shown in Figures ||-||. We compare the Pade protocol with results using the N/D protocol 
for three different values of the mass parameter M . 

Note that, for values of M below around 10 TeV there are no resonances at all in the 
N/D scenario. This is in accord with expectations based on Figures [2| to ||. Also, if 
M becomes too large then it leads to unusual behaviour of the amplitudes due to the 
dominance of the g{s) term which suppresses the amplitudes away from the region of 
resonances and can produce zeros in the individual partial wave amplitudes. The tail in 
the dotted line shown in Figure ^ is a consequence of such behaviour. Just discernable 
in Figure || is an isospin 2 scalar resonance just below 1.5 TeV in the N/D dotted curve. 



11 



500 



1000 1500 
M w (GeV) 



2000 



2500 



Figure 6: Parton level cross-section for Scenario B. We compare the Pade result with the 
N/D results as in Figure ||. 
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Figure 7: Parton level cross-section for Scenario C. We compare the Pade result with the 
N/D results as in Figure ||. 
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Figure 8: Parton level cross-section for Scenario D. We compare the Pade result with the 
N/D results as in Figure B. 
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Figure 9: Parton level cross-section for Scenario E. We compare the Pade result with the 
N/D results as in Figure ||. 
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Figure 10: Typical diagrams for signal and background processes: (a) signal; (b) VF+jets; 
(c) it 



5 Monte Carlo Simulations 



We have modified the Pythia Monte Carlo generator [15| to include the EWChL ap- 
proach using both Pade and N/D protocols. Signal samples containing the W^W^ final 
state (including all charge combinations) are generated using Pythia 6.146 with the 
Pade unitarisation scheme^. As a cross check, a sample with a 1 TeV Higgs was also 
generated using Herwig 6.1 j27j) . 



The dominant backgrounds are QCD tt production and radiative T^+jets, as illustrated 



in Figure |10|. These processes are implemented in the Pythia 6.146 and Herwig 6.1 
generators. To improve generation efficiency the minimum pt of the hard scatter is set to 
250 GeV for the W+ jets sample and to 300 GeV for the it sample [||. In addition to the 
hard subprocesses, the effects of the "underlying event" are simulated in both signal and 



background. Our default model in Pythia |28|| is obtained by setting a fixed minimum 



Pt cut off of p^ m = 3 GeV for secondary scatters. The default energy dependence of 
3 The code is available from the authors on request. 
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this cut-off has been explicitly turned off. No pile-up from multiple pp interactions is 
included. Other models, in both Herwig and Pythia, are discussed in Section ||, along 
with their effects. The leading order cross-sections are used to obtain rates and there is 
therefore a rather large degree of uncertainty, particularly in ti production, which is a 
pure QCD, dominantly gluon induced, process. NLO calculations [|^] suggest K-factors 
of order two are appropriate; the final word would come from measurements at the LHC 
itself. 



6 Extracting the Signal 



To identify semileptonic decays, we select first on the leptonically decaying W (elec- 
tron/muon and missing transverse energy), then on the hadronically decaying W (jet 
invariant mass, rapidity and transverse energy) and finally on the event environment 
(tagging jets at high rapidities, vetoing on central minijet activity). In all cases we have 
used only particles within a rapidity region of \rj\ < 4.5 to approximate the acceptance 
of a general purpose detector at the LHC. For clarity, we show just one signal sample as 
an example. The 1 TeV scalar resonance (scenario A) is chosen, since this has the lowest 
average Mww and therefore has a shape closest to that of the backgrounds. The other 
scenarios, while in general very like this sample, have a harder spectrum in the transverse 
momentum variables. The analysis follows the 1 TeV Higgs study of || quite closely for 
many cuts. However, we differ in the identification of hadronically decaying W bosons 
via the subjet method, in the top quark veto, in the cut on the transverse momentum of 
the hard system, and in details of other cuts; all of which are described below. 



6.1 Leptonic Variables 



Figure |ll| shows (a) the transverse momentum and (b) rapidity of the highest transverse 
momentum charged lepton for signal and background processes. The W+jets background 
is very similar to the signal in these distributions. Leptons from the it background are 
slightly softer and more central. Figure |ll|(c) shows the missing transverse momentum. 
Again, the it background is slightly softer than the other two samples. 

All leptons in an event are then combined one-by-one to give, if possible, a reconstructed 
W boson (to within a twofold ambiguity due to the unknown z component of the neutrino 
momentum). The transverse momentum of all these W candidates is shown in Figure 
|ll|(d). The signal has a harder distribution than both backgrounds. A selection cut is 
applied at 320 GeV on this distribution and in the case that more than one candidate is 
present, that with the highest transverse momentum is used. 



6.2 The Hadronic W Decay 

Figure |l~2l(a) shows the transverse momentum and (b) the pseudorapidity (ry) of the 
highest transverse momentum jet in the remaining signal and background samples. Jet 
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Figure 11: Leptonic variables for signals and background (a) transverse momentum of the 
highest pt charged lepton (e or fj,), (b) pseudorapidity of the same lepton, (c) 
missing transverse momentum and (d) the pt of the W candidate constructed 
from the lepton and the assumed neutrino. The area under the histograms 
is set to one to allow comparison of the shapes. A trigger cut at 80 GeV in 
the pt of the highest px jet and at 40 GeV in the highest px charge lepton is 
applied before making the plots, as well as a realistic rapidity acceptance. 
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finding is performed with the inclusive kx algorithm ijlj], and the E recombination scheme 
is used throughout. To reconstruct the W mass, the highest transverse momentum jet 
within the region \rj\ < 4 is selected. In the E recombination scheme the candidate W 



mass, Mj is then the invariant mass of this jet. Figure |T2](c) shows this distribution, 
with W mass peaks visible in the signal and in the tt sample, and a top mass peak also 
visible in the tt events. Cuts are applied at px > 320 GeV, and 70 GeV< Mj < 90 GeV. 
The results (after this cut and the leptonic cuts) are shown in the second and third rows 
of Table |. 

The jet is next forced to decompose into two subjets. The possibility of using subjets to 
reconstruct massive particles decaying to hadrons has been discussed previously [3"0|| . In 
this analysis we develop a new technique. The extra pieces of information gained from 
the subjet decomposition are the y cut at which the subjets are defined and the four- 
vectors of the subjets. For a genuine W decay the expectation is that the scale at which 
the jet is resolved into subjets (i.e. yp T ) will be 0(My/). The distribution of logipry/y) 
is shown in Figure |l2|(d). The scale of the splitting is indeed high in the signal and softer 
in the W+ jets background, where the hadronic W is in general a QCD jet rather than 
a genuine second W . A cut is applied at 1.6 < log(pxy/y) < 2.0. The effect of this cut is 
shown in the fourth row of the table. Whilst this is a powerful cut for reducing the W+ 
jets background, the effect on the tt background, which more often contains two real W 
bosons, is less marked. 



6.3 The Hadronic Environment 

To further reduce backgrounds, cuts must be applied to characteristics of the event other 
than those directly related to the decaying W bosons. 

Top quark veto In the remaining tt events containing a genuine leptonic W, the W 
will combine with a jet other than the hadronic W candidate to give a mass close to the 
top mass. This mass distribution for the leptonic W candidate combined separately with 
each such jet in the event is shown in Figure [0|(a). The top peak is clearly visible in the 
it sample. Any event with a mass in the region 130 GeV < M w j < 240 GeV is rejected. 
A similar distribution (not shown) is obtained by combining the hadronic W candidate 
with other jets in the event, and the same cut is applied. In combination these cuts are 
decribed as a "top quark veto", and their effect is shown row five of Table §. 



Tag jets In the WW scattering process the bosons are radiated from quarks in the 
initial state (see Figure 0(a)). The quark from which the boson is radiated will give 
a jet at high rapidity (i.e. close to the direction of the hadron from which it emerged). 
These jets are not in general present in the background processes and demanding their 
presence is therefore a powerful tag of the signal [13|. In this analysis we define a "tag jet" 
as follows. The event is divided into three regions of rapidity: "forward", i.e. forward of 
the most forward W; "backward", i.e. backward of the most backward W; and "central", 
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Figure 12: Kinematic variables for the hadronically decaying W candidate, (a) pr, (b) 77, 
(c) Invariant mass (d) PTy/y ■ The area under the histograms is set to unity 
to allow comparison of the shapes. 
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i.e. the remaining region, which includes both W candidates. A forward (backward) tag 
jet is defined as the highest transverse energy jet in the forward (backward) region. In 
Figure 13(b) the rapidity distribution of the tag jets with pr > 20 GeV is shown. Signal 
events display an enhancement at high \rj\ and a suppression at low \r]\, in dramatic 
contrast to the background processes, where most jets are central. For an event to be 
retained it must have a tag jet in both the forward and backward regions satisfying 
Pt > 20 GeV, E > 300 GeV and 4.5 > \rj\ > 2. The result of imposing this cut is shown 
in row six of Table |2|. The background is reduced by a factor of around fifty, at the cost 
of a loss less than two thirds of the signal. 



Hard px Figure |l3|(c) shows the px distribution for the "hard scattering" system com- 
prising the two tags jets and the two W candidates. For events surviving the cuts so far, 
the background events have a harder spectrum than the signal, since in the signal events 
this system is the complete result of a scattering between colinear partons, whereas in 
the backgrounds extra jets from hard QCD radiation may be picked up and/or missed. 
An upper cut is applied at 50 GeV, and the results are shown in row seven of Table |. 



Minijet Veto Finally, a cut which has been employed before in similar analyses ||, 11 



12|| exploits the fact that for signal events no colour is exchanged between the quarks 



which radiate the W bosons and the jets which are produced by the hadronically decaying 
W. This leads to a suppression of QCD radiation in the central region in the signal 
with respect to the background. However, significant activity is expected in all classes 
of event due to remnant-remnant interactions ("underlying event"). This activity can 
produce additional (mini)jets, and so it is important to choose a cut on additional jet 
activity which is robust against the large uncertainties in current understanding of the 
underlying event at the LHC. In this analysis minijets are defined as all jets apart from 
the hadronic W candidate with \r]\ < 2. Events are vetoed if the number of minijets with 
Pt > 15 GeV is greater than one. The distribution of the number of jets satisfying these 
demands is shown in Figure [O^(d). The result of applying this cut is shown in row eight 
of Table f| This cut is discussed further in Section || 



7 Analysing the signal 

7.1 Efficiency and Event Numbers 

Having applied the cuts described in the previous section, the WW mass distribution 
obtained is shown in Figure [F|(a) and (b) for all five signal samples discussed above. 
The resolution obtained in this variable is around 10 GeV, before any detector smearing. 
The efficiency is shown as a function of the true WW mass in (e). It rises from zero to 
6% between 500 GeV and 1.5 TeV, and is flat above this value. This efficiency includes 
the branching ratio for semileptonic W decays of around 15% . Excluding the branching 
ratio, the efficiency is around 40%. 
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Figure 13: (a) The mass distribution for the leptonic W candidate combined seperately 
with all other jets in the event other than the hadronic W candidate, (b) The 
rapidity distribution for tag jets (see text), (c) The transverse momentum 
distribution for the WW+ tag jets system, (d) The number of minijets (see 
text). 
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Figure 14: (a,b) Distribution of the reconstructed WW mass for signals and backgrounds 
separately. (c,d) Distribution of |cos0*|, the absolute value of the cosine of 
the centre-of-mass scattering angle for signals and backgrounds seperately. (e) 
Efficiency for signal events as a function of the true My/w and (f) | cos#*|. 
The errors reflect the statistics which would be obtained after approximately 
one year of running at the LHC, i.e. 100 fb _1 . 
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A key variable for distinguishing between scalar and vector resonances is the angular 
distribution of the scattered W pair in the WW centre-of-mass system. In Figure |l4](c) 
and (d) the distribution of | cos#*| is shown, where 6* is the angle between the scattered 
W and the incoming W direction, in the WW centre-of-mass frame. In (f) the efficiency 
is shown. The efficiency is very dependent on the mass distribution, since for the same 
transverse momentum, high scattering angles have high mass. This means that the 
transverse momentum cuts bias this distribution. However, this bias is well understood 
and could be corrected for in a final measurement using a two-dimensional correction in 
mass and angle regardless of the input distribution. 



7.2 Simulated Measurement 

If it is assumed that the backgrounds can be well constrained from developments in the- 
ory, measurements at the Tevatron and HERA over the next few years, and measurements 
at the LHC in other kinematic regions, then the statistical error on an extraction of the 
My/w and | cos#*| distributions can be estimated by adding the statistical errors on the 
signal and background distributions in quadrature. Under this assumption, a simulation 
of an expected measurement of the differential cross-section da / dMy/w after 100 fb _1 
of LHC luminosity is shown in Figure |l5|(a), (c) and (e). The scenarios containing reso- 
nances are distinguishable above the background, and are also distinguishable from each 
other due to their different resonant masses. In Figure [l^(c) the double resonance sample 
(D) is shown, with two peaks clearly measured. Also shown (in all three figures) is the 
continuum model (E). 

The expected measurement of the differential cross-section da/d\ cos #*| after 100 fb _1 



of LHC luminosity is shown in Figure |l5|(b),(d) and (f) for My/w > 750 GeV. The inter- 
mediate mass vector and scalar resonances have the expected behaviour, with the vector 
rising towards high | cos 6** | and the scalar being flat. In Figure (l^(d) the distribution 
for the double resonance model is shown in two mass bins: 750 < My/w < 1200 GeV 
and My/w > 1200 GeV. With the high statistics generated (corresponding to a very high 
integrated luminosity), the lower mass resonance can be seen to be a scalar whilst the 
higher mass is a vector. However, within the simulated errors the measurement of the 
spin of the lower mass resonance would be marginal. 



8 The Underlying Event 

One of the more uncertain aspects of the analysis is the understanding of the so-called 
"underlying event". This is defined here as particle and energy flow in the event associated 
with the same proton-proton interaction but incoherent with the W production process. 
Hence we explicitly exclude from our definition the effects of multiple pp interactions 
in the same bunch crossing, any detector effects such as those associated with noise or 
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Figure 15: Measurement expectation after 100 fb of LHC luminosity at 14 TeV cm 
energy. (a,c,e) da / dMww and (b,d,f) da/d\ cos0*|. (d) shows da/d\ cosf?*| for 
the high and low mass subsamples for the double resonance model, separated 
by a cut at 1200 GeV. 
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Cuts 



Generated 



p T (Lep. W )> 320 GeV 
and 

p T (Had. W) > 320 GeV 



70 GeV < M(Had. W 
< 90 GeV 



1.6 < log(;p T x ^ ) < 2.0 



Top quark veto 
(see text) 



Tag jets 
p T > 20 GeV, E > 300 GeV 
(see text) 



Hard pt < 50 GeV 



Minijet veto 
Pt > 15 GeV, see text 



Efficiency 



A:100% 
B:100% 
C:100% 
D:100% 
E:100% 



A:ll% 
B:ll% 
C:10% 
D:10% 
E:10% 



A:6.7% 
B:6.2% 
C:5.8% 
D:5.6% 
E:5.8% 



A:4.7% 
B:4.4% 
C:4.1% 
D:4.0% 
E:4.1% 



A:4.3% 
B:4.0% 
C:3.8% 
D:3.6% 
E:3.8% 



A:1.6% 
B:1.5% 
C:1.4% 
D:1.3% 
E:1.4% 



A:1.5% 
B:1.5% 
C:1.4% 
D:1.3% 
E:1.4% 



A:1.5% 
B:1.5% 
C:1.4% 
D:1.3% 
E:1.4% 



Signal 
a(fb) 



72 
104 
44 
113 
47 



8.2 
11 
4.4 
11 
4.7 



4.8 
6.4 
2.6 
6.3 
2.7 



3.4 
4.5 
1.8 
4.5 
1.9 



3.1 
4.2 
1.7 
4.1 
1.8 



1.1 
1.6 

0.63 
1.5 

0.67 



1.1 
1.5 

0.61 
1.4 

0.65 



1.1 
1.5 

0.61 
1.4 

0.65 



tt 
a (lb) 



W+Jets 

<7(fb) 



Pythia 
18,000 | 65,000 

Herwig 
14,000 53,000 



Pythia 
910 | 4400 
Herwig 



750 



3600 



Pythia 
56 | 700 
Herwig 



52 



480 



Pythia 
28 | 78 
Herwig 



27 



66 



Pythia 
3.2 j 52 
Herwig 



3.4 



43 



Pythia 
0.030 | 0.38 
Herwig 



0.082 



0.42 



Pythia 
0.020 | 0.32 
Herwig 



0.048 



0.37 



Pythia 
0.013 | 0.24 

Herwig 
0.048 0.36 



Sig/B 



8.7 x 10" 4 
1.3 x 10~ 3 

5.3 x 10~ 4 

1.4 x 10~ 3 
5.0 x 10~ 4 



1.5 x 10" 3 
2.1 x 10~ 3 
8.3 x 10~ 4 
2.1 x 10~ 3 
8.8 x 10~ 4 



6.3 x 10" 3 

8.4 x 10~ 3 
3.4 x 10~ 3 
8.3 x 10~ 3 
3.6 x 10~ 3 



3.2 x 10~ 2 

4.3 x 10~ 2 

1.7 x 10~ 2 
4.3 x 10~ 2 

1.8 x 10~ 2 



5.6 x 10" 2 
7.5 x 10~ 2 
3.0 x 10~ 2 
7.3 x 10~ 2 
3.2 x 10~ 2 



2.7 
3.8 
1.5 
3.6 
1.6 



3.2 
4.5 
1.8 
4.3 
1.9 



4.3 
6.0 
2.4 
5.6 
2.6 



Table 2: The effect of cuts on the signal and background samples. A: 1 TeV scalar, B: 
1.4 TeV Vector, C: 2 TeV Vector, D: Double Resonance and E: Continuum. 
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pile-up, and hard QCD radiation associated directly with the hard scatter. The first 
two of these are not simulated here and controlling and understanding them requires 
detailed experimental work. The third is simulated to leading-logarithmic accuracy in 
both Pythia and Herwig. While this simulation should and probably will be improved 
in the future, for now it is considered adequate. 

The remaining activity can be characterised as interactions between the proton remnant 
systems. It is important because it is largely independent of the hard scattering process, 
and therefore contributes to minijet activity in both signal and background, degrading 
the effectiveness of the minijet veto. In addition, underlying event activity contributes to 
the observed W width and the position of the mass peaks in a highly model-dependent 
way. 

In Figure |l6](a) and (b) the jet mass distribution and the iog(pry/y) are shown again 



(as in Figure 12(c) and (d)) for the signal events (1 TeV resonance, sample A) using 
our default underlying event model. In addition, several other underlying event models 
are shown. In PYTHIA, we turn off multiparton interactions (sample Al), and turn on 
the default model (sample A2) which has a p™ n of 2.89 GeV at LHC energies. Also 
shown are three samples of 1 TeV Higgs events generated using HERWIG. These have 
no underlying event (sample A3), soft underlying event (sample A4) and multiparton 
interactions generated with fixed p™ zn = 3.0 GeV (sample A5)fl3l]. The width of the W 
mass peak is much greater in general for those samples which include an underlying event. 
Whilst the PYTHIA multiparton interaction models and the HERWIG soft underlying event 
are fairly consistent with each other, the Herwig multiparton interaction model gives a 
very different distribution. However, the log(pTy/y) is similar for all models, implying 
that this cut should be robust against such uncertainties. 

For the same samples, the minijet px distribution and the number of minjets passing 
the 15 GeV cut, which we introduced in the analysis of Section [6^, are shown in Figure 
|l6|(c) and (d), with absolute normalisation. In contrast to the W mass distribution, in 
these distributions the Herwig multiparton interaction model is close to the PYTHIA 
multiparton models, whereas the soft underlying event model is closer to the models 
without underlying event. The pr distribution is very steeply falling, and is sensitive 
to the underlying event below around 20 GeV. Thus, there is sensitivity in the number 
of jets at 15 GeV, and this would become worse for lower choices of cut. Lowering the 
cut further without introducing large uncertainties requires a better knowledge of the 
underlying event than is currently available. 

If the no underlying event model is used in PYTHIA (sample Al), the signal/background 
for the scenario A is 8.0. However, for all other cases (models A, A2-A5) the ratio 
is between 2.5 and 4.0. Data from the Tevatron and photoproduction at HERA (see 
for example [^] and references therein), strongly disfavour models without underlying 
event (Al, A3) and are generally more consistent with the other models considered 
here (though none provides a perfect description). However, further work is needed on 
constraining these models to improve confidence in the extrapolation to the LHC. At 
present a systematic error of 40-50% would have to be assigned to the measurement from 
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this source alone. 



9 Summary and Conclusions 

A major goal of the LHC is to extract the WW — > WW cross-section as accurately 
as possible to the highest centre-of-mass energies in order to shed light on the nature 
of electroweak symmetry breaking. We have performed a study of the WW — ► WW 
scattering cross-section in the scenario that there is no new physics below the TeV scale 
using the formalism of the Electroweak Chiral Lagrangian extended by the imposition 
of unitarity constraints. Two different unitarisation protocols are used: Pade and N/D. 
These protocols determine the behaviour of the scattering cross-section into the TeV 
regime and they typically predict the emergence of new vector and/or scalar resonances. 
We have performed a detailed comparison of these two unitarisation methods. 

We have implemented the physics of the unitarised Electroweak Chiral Lagrangian in 
a realistic general-purpose Monte Carlo (Pythia). The semi-leptonic decay mode of 
the final state W pair has been studied at the final state particle level with detector 
acceptance cuts but no smearing. We have considered five different physics scenarios 
which are representative of the different types of physics which we might reasonably 
expect at the LHC. The principal backgrounds come from W+ jet and tt production, 
and we consider these backgrounds using both the PYTHIA and HERWIG Monte Carlos. 
A new method for identifying hadronically decaying W bosons is introduced which we 
expect to be useful more generally in the identification of hadronically decaying massive 
particles which have energy large compared to their mass. Other new features include 
a top quark veto and a cut on the transverse momentum of the hard subsystem. In 
addition, the established tag jet and minijet veto cuts are applied. The results are cross- 
checked with Herwig using a simulation of a 1 TeV Higgs boson for the signal. The 
effect of uncertainties in the underlying event leads to a model dependent systematic 
error of 40-50%. New data from Tevatron and HERA should help to reduce this before 
the LHC turns on. 

The results compare very well with previous Higgs search studies in the semi-leptonic 
channel. Over a wide range of parameter space signal/background ratios of greater 
than unity can be obtained, and the cross-section can be measured differentially in the 
WW centre-of-mass energy within one year of high luminosity LHC running (100 ft) -1 ). 
Vector and scalar resonances up to around 1.5 TeV may well be observable, and their 
spins measureable. Detailing the exact regions of sensitivity, as well as verifying the 
improvements in signal/background arising from the new cuts, requires a more detailed 
simulation of the LHC general purpose detectors. 
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Figure 16: The effect of the underlying event, (a) Hadronic W mass, (b) pr^/y , (c) 
Number of minijets and (d) the pr distribution of the minijets. 
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